Quantum control and manipulation of donor electrons in Si-based 

quantum computing 

Maria Jose Calderon 
Instituto de Ciencia de Materiales de Madrid (CSIC), 
Cantoblanco, 28049 Madrid (Spain) 

Andre Saraiva and Belita Koiller 
Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Caixa Postal 68528, 21941-972 Rio de Janeiro, Brazil 

Sankar Das Sarma 
Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, MD 20742-4111 
(Dated: July 21, 2009) 

Abstract 

Doped Si is a promising candidate for quantum computing due to its scalability properties, 
long spin coherence times, and the astonishing progress on Si technology and miniaturization in 
the last few decades. This proposal for a quantum computer ultimately relies on the quantum 
control of electrons bound to donors near a Si/barrier (e.g. Si02) interface. We address here 
several important issues and define critical parameters that establish the conditions that allow the 
manipulation of donor electrons in Si by means of external electric and magnetic fields. 
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I. INTRODUCTION 



In the last few decades, the density of transistors in a chip has been consistently duplicated 
every two years, as given by Moore's law. This demand is pushing the semiconductor 
devices fabrication (in particular, for Si) towards the atomic limit in which the number of 
dopants is so small that their exact position and distribution affects the performance of 
a transistor.^ There is therefore a great interest in achieving an atomic control over the 
position of impurities in Si.^'^ This interest has been reinforced by the proposal of using 
shallow donors (typically P or As) as spin qubits for a Si-based quantum computer,-!^ and 
related ones that followed,- where qubit operations are performed by manipulating the donor 
electron with an external electric field provided by local gates. Spin is a natural candidate 
for a qubit, particularly in Si where very long spin coherence times have been measured (at 
least ~ 1 ms in natural Si).-i^>^i^ Spin coherence time has to be orders of magnitude longer 
than the time required for operations, and enough to allow for quantum error correction.— 

In the doped Si-based quantum computer,-*^ schematically shown in Fig. [T], single-qubit 
operations (rotations of the spin state on the Bloch sphere) are accomplished with an AC 
magnetic field in resonance with the level splitting of the hyperfine coupled nuclear spin- 
electron spin system. The hyperfine coupling is proportional to the probability density of 
the electron wave-function on the donor. Therefore, the on and off resonance condition 
is controlled by applying electric fields to manipulate the electron wave-function. Two- 
qubit operations could in principle be performed by letting electrons in neighboring donors 
interact via a transient exchange coupling during a specified period of time.^^ Combinations 
of these exchange mediated operations and single spin rotations configures a controlled-NOT 
(CNOT) gate.^^ With the one qubit rotations and the CNOT gate any operation on n qubits 
may be performed.— 

We have considered the problem of a donor a distance d from a Si/Si02 (001) interface and 
studied the manipulation of the donor electron by electric and magnetic fields. Our study 
applies for cases when d is large enough so that, under an applied field, abrupt ionization (via 
tunneling) takes place.— When no external fields are applied, the only attractive potential 
felt by the electron is the donor Coulomb potential and, at the low temperatures relevant 
here, the electron remains bound to the donor [see Fig. [2ta)]. Under an electric field {F) 
applied perpendicular to the interface, a triangular potential well is formed at the interface. 
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FIG. 1: Basic scheme of the doped Si-based proposal for quantum computing. ^ It consists of an 
array of donors in Si separated by a distance R and located at a distance d from a Si02 barrier. The 
donor electrons are manipulated by applying electric fields by local gates. The A gates control the 
relative position of the electron with respect to the donor affecting the hyperfine interaction between 
their spins (allowing 1-qubit operations). The J gates control the overlap between neighboring 
electrons, allowing to perform the exchange pulse involved in 2-qubit operations. 



as shown in Fig. Mjo). For fields larger than a certain minimum value, bound states exist at 
the triangular well. When F reaches a characteristic strength F^. such that the bound state 
at the interface is degenerate in energy with the bound state at the donor, tunelling between 
the two wells (the donor Coulomb potential, and the triangular well at the interface) is 
possible. In this way, we can reversibly manipulate the position of the electron between the 
donor and the interface by oscillating the electric field around F^ [namely, going back and 
forth from situation in Fig. [2](b) to situation in Fig. EJ^c)]. Note that the electron, when 
at the interface, may still remain bound to the donor along the xy-plane, and hence, the 
electron may be confined in all three space directions upon certain conditions that will be 
described below. We simplify the problem by assuming that the Si02 barrier is infinite. 

The outline of this paper is the following. In Sec. [TTl we describe our model and discuss 
the peculiarities of the conduction band of Si, that has a multivalley structure. We also 
describe our basis set consisting on the lowest uncoupled donor and interface states. In 
particular, we point out the conditions that guarantee the confinement in all 3 directions 
of space of the electron at the interface. In Sec. IIIII we discuss the conditions that allow 
the manipulation of the donor electron to be possible, in particular:— (i) we calculate the 
value of the characteristic electric field at which shuttling of the electron between donor and 
interface may occur; and(ii) we estimate the tunneling time of this process. We also briefly 
discuss the effect of a magnetic field applied parallel to the electric field.— In Sec. HV] we 
discuss how the 2-qubit operations could be performed with the electrons at the interface 
stated and show new results for the calculation of the exchange calculated with an improved 
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FIG. 2: (a) If no electric field is applied the electron only sees the donor potential and is bound 
to it. (b) When an electric field F, perpendicular to the interface, is applied, a triangular shaped 
potential forms at the interface. For a certain value of F a bound state exists at the interface, 
(c) When F ^ the states at the interface and at the donor are degenerate and tunneling 
between donor and interface may occur, (d) Donor potential at the interface along the xy-plane 
V{x) = • This attractive potential causes the electron not to spread along the interface 

when it is ionized from the donor. 



Heitler-London method. We finish with a discussion and conclusions in Sec. |Vl 



II. MODEL 

The wavefunction of a donor electron in a semiconductor can be written as an expansion in 
terms of the Bloch waves close to the bottom of the conduction band.— The conduction band 
of Si has six minima (valleys) located in the (100) directions at a distance kg = 2tt 0.85/ asi 
{asi = 5.4 A the lattice parameter of Si) from the T point. Therefore, combinations of Bloch 
waves from the 6 valleys have to be considered, and the donor ground state wave-function 
would be 6-fold degenerate. This degeneracy is lifted when the coupling between Bloch 
waves from different valleys due to the singular donor coulomb potential is taken into account 
(the so-called valley-orbit coupling). The resulting ground state is non-degenerate and it 
is a symmetric combination of the six valleys (symmetry Ai).^*^ Within the effective mass 
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approximation, only the Bloch functions at the positions of the conduction band minima 
are involved, and the ground state of the electron at the donor is written^ 

n = ^J2^Di^^M^,rD), (1) 



where F^{r) are envelope functions (/i = +x,-x,+y,-y,+z,-z), and 0^(r, r^) = 
^Bioch^"*''^ ""^ ~ Uf^{r)e'^^^'''--^~'^°\ Here, rn is a reference site, which represents a pinning 
site for the Bloch waves in the superposition state. In the case of a single Bloch state, td 
contributes with an irrelevant phase, but for superposition states it leads to interference 
effects. For an isolated impurity, r^j is naturally chosen to be the position of the donor 
(td = 0). 

For the description of the interface state it has to be taken into account that the (001) 
interface breaks the degeneracy of the 6 valleys, making the x, -x, y, and -y valleys much 
higher in energy than the z and -z valleys.— (A tensile strain has the same effect on the 
valley degeneracy).— Therefore, the ground state at the interface only includes the Bloch 
states from the z and -z valleys. We write it 

vl/f = -^Fiir) [0,(r, zi + d)± 0„,(r, Zi + d)] (2) 

where the envelope function Fj is obtained variationally assuming the form: 

Fj = ^(z + df e-^^^+'^^Z^ X 4= e-^"'"'^ (3) 
2v 3 V TT 

with a and (3 taken as variational parameters, which are the same as for the single valley 
approximation.— The parameter 1/a is related to the width of the wave-function at the 
interface in the z-direction while 1/(3 gives the width of the wavefunction along the xy 
plane.— >^ For the isolated interface, we assume the pinning site zi is exactly at the interface. 

The value of depends on how far the donor is from the interface. When the electron 
is at the interface it still feels the attractive potential of both the donor and its image [see 
Fig. [2] (d)]. Typical values of range from 6 nm for d = Q nm to 22 nm for (i = 40 
nm. If we want to be able to shuttle the electron between the interface and the donor 
reversibly, we need the electron to be confined at the interface, and not to spread forming a 
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two-dimensional gas. This gives us a limit for the maximum planar density of around 10^° 
cm-2.i^^ 

The fact that the 6-valley degeneracy is broken at the interface plus the fact that the 
smooth potential between the donor and the interface does not couple valleys in different 
directions, allows us, in a first approach to the problem, to consider only the z and -z valleys 
for the state at the donor. This simplification would be exact for strained Si. Therefore, it 
makes sense to use 

vl/g = -^Fnir) [0,(r, zn) ± 0-.(r, Zn)] (4) 

with 

Fn(x{z + d)e-V^^^^+^ , (5) 

the donor envelope function centered at r = 0. This is a hydrogen-like envelope with an 
anisotropic shape that arises due to the conduction band effective mass anisotropy in Si, 
multiplied by the factor {z + d) that guarantees that the wavefunction fulfills the boundary 
condition for an infinite barrier at ^ = —d. We are interested here in the situation where the 
donor is relatively far from the interface, i.e, the ground state wavefunctions in each well do 
not overlap significantly, and in this range the variationally determined Bohr radii a and b 
are the same as for a donor in bulk (a = 2.365 nm and b = 1.36 nm).'^^ 

Note that, for each isolated well, the value of the splitting between the symmetric and 
antisymmetric combinations of the Bloch waves and which one is the ground state depend 
on the value of the pinning positions z^ or zj, as illustrated in Fig. [51 

Let us write now the full Hamiltonian for a donor (at 2; = 0) in Si a distance d from a 
(001) interface with Si02. The boundary problem between the two semi-infinite dielectrics 
is considered by including the image charges for the donor and for the electron (see Fig. H]). 
A uniform electric field, perpendicular to the interface, is also included. The resulting 
Hamiltonian is,— in rescaled atomic units a* = h'^€si/m±e'^ = 3.157 nm and Ry* = 
m±e^/2ffe% = 19.98 meV, 

H = T + KeFz - ^ + V^^'""^" + V;^"^^se + ^ (6) 

where T = — ^ — ^ — 7^, 7 = 'm'±/'m\\, "mi = 0.191m, my = 0.916, k = 3.89 x 
10^''e|j (m/m_|_)^ cm/kV, and the electric field F is given in kV/cm. The next terms are the 
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FIG. 3: Typical plots of the shapes of the envelope functions at the donor (a) given by Fq Eq. ([5]), 
and at the interface (b) given by Fj Eq. ([3]). -F/ is much more extended that F^, in the xy-plane 
~ 10 — 20 nm, depending on d, — ill while the Bohr radii a = 2.365 nm). 
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FIG. 4: We consider Si and Si02 to be two semi-infinite layers. We include the image charges D' 
and e'. For our particular value oi Q = —0.5 we see that the sign of the images is the same as that 
of the original charges. Therefore, the donor image potential is attractive and the electron image 
potential is repulsive [see the expressions for V^™'^^'' Eq. ([7]) and y^^^^^ Eq. dS])]. 
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FIG. 5: The expectation values for the energies corresponding to the symmetric (sohd hne) and 
antisymmetric (dashed Une) combinations of the z and -z Bloch waves depend on the value of the 
pinning site Zj, varying as cos(2/co-2j). For Zj = the ground state is the symmetric combination 
and the splitting between the two combinations is maximum and equal to 2V^^. But for very small 
displacements of the pinning site (of less that 1 A) the splitting changes enormously. The dotted 
horizontal line corresponds to the ground state energy within the single valley approximation ! ^^i^'' 

Q is a function of the dielectric constants Q = (esi02 ~esi)/(esi02 +fsi)) with esi = 11.4 
and esi02 = 3.8. With this Q we get that the donor image potential is attractive while the 
electron image potential is repulsive. 

The last term in the Hamiltonian Eq. ([6]) describes valley-orbit effects, namely the cou- 
pling between different valleys due to the singular nature of both the donor (D) and the 
interface (/) potentials. These couplings are quantified by the parameters Vj^ and Vj^. Vj^ 
is known from the splitting of the IS manyfold of the isolated donor spectrum = —1.5 
meV. We use this value all throughout. Vj^ is not precisely known, although estimated to be 
in the order of 1 meV,— and probably dependent on the interface quality.— We consider 
it as a parameter ranging from to —10 meV. The results described below for a finite value 
of the interface valley orbit coupling do not depend qualitatively on the particular value or 
sign of Vi for \ Vl\ > 0.02 meV. 

In our two-valley formalism, assuming that d is large enough so that no strong hy- 
bridization occurs between the donor and interface states, the problem may be restricted to 
the basis set of the lowest uncoupled donor and interface states, namely {-^D^Bioch^"*'^''^^' 
^i'^Bitch^''"'^''' ^/^Bioch^"*'""'', ^/^moch^''""''}' leading to the Hamihonian matrix 
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where Ed = {Ed\H\Ed) , Ej = {Fj\H\Ej), and Hjd = {Fj\H\Fd) are the same as the 
single-valley matrix elements.— 

III. MANIPULATION OF A SINGLE DONOR ELECTRON CLOSE TO A SI/SIO2 
INTERFACE 

We solve the Hamiltonian Eq. for different (i's and a range of electric fields F. At each 
point, the pinning sites zj and zd are determined variationally such that the total ground 
state energy is minimized. For F << F^, the electron is at the donor and the pinning is 
on the donor, to be precise, abs[cos(2kQZD)] = 1 while a6s[cos(2A;o(-2/ + d))] < 1. In the 
opposite limit, F » Fc, the electron is confined at the interface and so is the pinning, 
namely, a6s[cos(2A;o(-2/ + d))] = 1 while abs[cos{2kQZD)] < 

The variational solution of the Hamiltonian Eq. leads to four eigenvalues that, as a 
function of the electric field, show certain level anticrossings. As an illustration, two typical 
cases are shown in Fig. O— When V^^ = 0, the two states at the interface are degenerate 
and only two anticrossings, involving three levels each, occur. For the more general case 
of Vj^ 7^ 0, there are typically four anticrossings, involving two levels each. For particular 
values of d, corresponding to cos{kQd) = 0, the coupling between levels is so strong that one 
anticrossing can involve all four levels. 

To determine the characteristic field Fc at which the 'shuttling' of the electron between 
donor and interface may occur and the times involved in the process, we need to look at 
the anticrossing between the lowest eigenvalues. The value of the electric field at which the 
anticrossing happens is Fc and the gap between levels at anticrossing gives an estimate for 
the tunneling time r ~ h/Egap. In Fig. [7] we show the value of Fc versus d for the single 
valley approximation, the 2- valley model with Vj^ = 0, and Vj^ = — 1 meV. The three curves 
show the same qualitative shape. The small shift between them arises due to the changes 
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produced in the ground state energies when the valley orbit coupling is included (this change 
is illustrated, with respect to the single valley energy, in Fig. [5]). Therefore, the multivalley 
structure of the conduction band of Si does not significantly affect the value of F^. This 
was, in fact, expected from previous tight-binding results that took into account the six 
valleys.— Very recent experiments on a Si FinFET with a single donor^^ show the same kind 
of behavior for versus d that we find here. 

Now we look at the value of the gap between levels calculated at anticrossing: The results 
are shown in Fig. [HI For the 2- valley model with degenerate states at the interface {V^^ = 0) 
the result is qualitatively similar to the smooth decrease already observed in the single valley 
approximation.— This is due to the fact that only one of the degenerate levels at the interface 
is coupled to the lowest level of the donor, leading to an effectively single valley result. 
However, a qualitative striking difference is observed when a finite valley orbit coupling is 
present at the interface: The gap versus d shows oscillations that go as abs[cos{kQd)], namely, 
the oscillations are not commensurate with the lattice. The tunneling times we get within 
the single valley approximation range from subpicosecond to nanoseconds depending on d 
(larger d correspond to longer vales of r). The times corresponding to the multivalley result 
can be much larger for d^s such that cos{kQd) ~ 0, but the statistical weight of such points 
(see the stars in Fig. [8]) is relatively small. Therefore, in an ion implanted sample, where the 
control on the donor depth positioning is limited, most donors would have tunneling times 
comparable to the single valley results, while others would take much longer. 

We have also analyzed the effect of a magnetic field applied parallel to the electric field.— >^ 
The magnetic field shrinks the electron wavefunction in the direction parallel to the interface. 
The effect is much stronger for less confined wavefunctions. Consequently, the wavefunction 
at the interface is much more affected than at the donor. The increase in confinement is 
concomitant with an increase in energy, which is hence significant at the interface while 
very small at the donor for magnetic fields on the order of a few Tesla. This magnetic field 
induced shift in energy has an interesting consequence: Starting from an electric field just 
above Fc (so the electron is at the interface), the application of a magnetic field can push 
the electron back to the donor.— The manipulation of electrons using simultaneously electric 
and magnetic fields can help us distinguish donor electron from other moving charges in real 
systems.— 
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FIG. 6: (a) Eigenvalues of the Hamiltonian in Eq. ([9]) for d = 5.99a* and Vj^ = 0. There are 
two anticrossings involving three levels each (the intermediate level remains uncoupled to the other 
two, which display the anticrossing) . (b) Eigenvalues of the problem for d = 5.94a* and V^^ = — 1 
meV. The four levels show anticrossings in pairs. The relevant anticrossing for the determination 
of Fc and the tunneling times (or, equivalently, gaps) is the one involving the ground state. The 
gaps at anticrossing are small at the energy scale of the figure. 



IV. EXCHANGE BETWEEN ELECTRONS AT THE INTERFACE 

We have discussed above how the multivalley nature of the conduction band of Si produces 
oscillations in the tunneling time involved in the shuttling of an electron between donor and 
interface. Similar oscillations have been found in the exchange interaction between electrons 
in neighboring donors as a function of the distance R between them.— 

Since the donor electron is much less confined at the interface than at the donor well, 
the transient exchange interaction between donors very distant from each other in the bulk, 
where they do not interact, could be performed by pulling (with an electric field) the two 
corresponding electrons towards the interface, where they would interact. Besides, if the 
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FIG. 7: Characteristic field at which tunnehng between the donor and the interface state may 
occur. The three curves (corresponding to single valley approximation, 2-valley with Vj^ = 0, and 
2-valley with Vj^ = — 1 meV) have the same qualitative behavior. 
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FIG. 8: Gap (related to tunneling times as r ~ h/Egap) versus d for single valley approximation, 
= 0, and = — 1 meV. The first two curves smoothly decay with d. However, the inclusion of 
a finite valley-orbit coupling at the interface leads to oscillations on the gap that go as abs[cos{kQd)\. 
The stars correspond to values of d given by monolayer steps of 1.35 A. 

two neighboring donors are at the same distance d from the interface, we do not expect to 
get oscillations on the exchange versus interdonor distance R. The 2D interface potential 
formed by two donors separated a distance R from one another and at the same depth d 
from the interface has the double- well shape, depicted in Fig.Oalong the axis (x) connecting 
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the well minima, as given by the expression 



V{x) 
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The distance R between donors has to be large enough so that the single electron wave- 
function is below the barrier between the wells, to guarantee that the electron shuttling 
between donor and interface is reversible. In Fig. [10] we show the exchange calculated for 
donor pairs at three different distances from the interface as a function of the interdonor 
distance. The calculations were performed in 2D, using the Heitler-London approximation 
with an improved hybrid variational wave-function which is gaussian in the center, where 
the parabolic approximation for the potential is valid and has an exponential decay at long 
distances of the well minimum, where the potential saturates and deviates significantly from 
the parabolic behavior,— 



This ansatz involves five parameters, three of which (Ai, y42, and rj) are obtained from the 
boundary conditions and normalization, and we are left with two variational parameters 
and /i). We note that r] is obtained from 77 = 2/5^/i. The adopted matched variational 
wavefunction gives an exchange coupling which fits well with the asymptotic values given 
by more rigorous calculations^ and which is orders of magnitude larger than the one we 
previously calculated using a pure gaussian variational form for the electron wavefunction 
at the interface.— The exchange we are getting is of the order of the one estimated for gated 
quantum dots in GaAs.— 

V. DISCUSSION AND CONCLUSION 

We have discussed here some of the relevant parameters that would allow elementary 
operations involved in the manipulation of donor electrons in Si close to a Si02 interface for 
quantum computing applications. For our calculations we have considered an ideal semi- 
infinite Si02 barrier. However, in real devices^'^'^ metallic gates are located on top of a 
relatively thin (few nm wide) Si02 barrier. It seems then that a more realistic description 
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FIG. 9: Donor pair potential at the interface, along the axis connecting the wells minima as given 
by Eq. pU|) for R = 28 nm and d = 6.3 nm. 
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FIG. 10: Exchange calculated within the Heitler-London approximation using a matched varia- 
tional wavefunction that correctly includes exponential decaying tails. 



of the devices should take into account the metallic gates. There exists a calculation similar 
to ours in the opposite limit of an infinitesimal Si02 barrier separating the Si layer from 
the metallic gate.— In that case, Q = 1 and the image charges have opposite sign to the 
real charges. In particular, the donor image potential will be repulsive instead of attractive 
and this can seriously deteriorate the confinement of the electrons along the xy-plane at 
the interface. For a realistic case of a Si02 barrier width of ~ 2 nm and an interface-donor 
distance (1 = 4 — 5 nm,— we would be in between the two extreme cases of a semi-infinite 
Si02 barrier considered here and an infinitesimal one considered in Ref. is^. In particular, 
the width of the wavefunction at the interface l/f5 would be larger than reported here. By 
comparing our 1/(3 with the results reported in Ref. Issi [see \/~(? ja^ in the inset of Fig. 4] we 
can estimate at most a doubling of 1//5 when metallic gates are included. This pushes the 
upper limit for the donor planar density towards smaller values, still keeping our qualitative 
conclusions valid. 



14 



Another relevant issue is related to the modification of the spin coherence time T2 when 
donors are located close to a Si surface. Different experiments^!^ have revealed that T2 
gets smaller than in bulk when donors are located close to an interface: bulk values reach 
jibuik ^ gg ^j-^jig close to a surface they range from T|™^ ~ 0.1 — 1 ms^ to T™"^^ ~ 1 /is^ 
for different samples. T2 depends on the quality of the interface (it is longer for a hydrogen 
passivated surface than for a Si/Si02 interface^) and the distance from the donor to the 
interface d. Some effects that may be responsible for the decrease of coherence times are the 
magnetic noise produced by impurities and other defects (like Pf, centers) at the interface,— 
and/or the recombination of donor electrons on the P;, centers.— These observations imply 
that increasing the interface quality may be crucial for the practical implementation of a 
Si-based quantum computer. 

In summary, we have described the basic conditions that would allow the manipulation 
of donor electrons in Si by external fields, taking into account the multivalley structure of 
the conduction band of Si. Experiments are getting close to the limit of isolated donors^^i^ 
where the operations discussed here are relevant. In fact, qualitative agreement for the 
behavior of Fc versus d has already been found in samples with isolated donors in a Si 
FinFET,— in a situation where the donor is closer to the interface than we considered in 
the present study. Both situations (small and large values of d as compared to the effective 
electronic confinement lengths) are important, and should be carefully investigated in the 
context of applications of donor electrons as carriers of quantum information through its 
charge or spin degrees of freedom. 
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